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- ABSTRACT 


A deep ocean capsule is proposed capable of sensing 
pressure, temperature, flow and direction for a study of 
tides and internal waves of the oceans. One of the major 
requirements for this capsule is the detection and recording 
of unusual events. To meet this need a very low power, 
highly reliable digital detector has been designed. COS/MOS 
logic with maximum use of MSI and LSI was used in the de- 
Sign. The detector was simulated uSing proposed word 
lengths and fixed point arithmetic which would be utilized 
in the ensuing hardware. The fixed point simulation was 
necessary to verify the problems in truncation and round-off 
errors and to provide necessary design criteria to ensure 


proper filter functioning and integral detection. 
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i. INERODUCTION 


An unusual event detector is to be a part of a one-year 
Geep ocean experiment [ Ref. 1]. The sensors anid veeera= 
ing system are all digital and have low power and high 
reliability as overriding design constraints. 

The detector is to have several distinct subsystems as 
“shown in Figure 1: a spurious data rejection gating system; 
a Semi-diurnal and diurnal tide rejection filter; anda 
saGuare and integrate fading memory detector. Upon determi- 
nation of an unusual event characterized by significant 
energy in the spectrum from 0.25 CPH to 30 CPH [Ref. l], 
various permanent memory and accoustic telemetry links may 
be initiated. The sensor signals are to be sampled at one- 
minute intervals and stored in a two-hour buffer memory. If 
no event is occurring, the signals are lost at the end of 
two hours. If an event is occurring, the buffer memory is 
tied anto the permanent memory. 

The sourious data rejection system is accomplished by 
tapping tne buffer memory at the half-hour points and taking 
first differences of adjacent signals. If the differences 
exceed a certain tolerance level, a spurious data point or 
glitch is assumed, and the first difference is rejected and 
is replaced with an appropriate value. These five first dif- 
ference .are then combined with appropriate weighting to 


achieve -) anti-tide filtering. The resulting signal is 
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then squared and fed into an integrator with a fading 
memory time constant of one-half hour. That is, the detec- 
tor will be predominantly sensitive to events that have 
accumulated over the past half-hour. 

Because the system is to operate unattended at the 
3,000-fathom level for one year, power and reliability are 
tantamount design considerations. For this reason, COS/MOS 
IC's are to be used with maximum utilization of existing LSI 
and MSI components. Based on these considerations, the 
following hardware design is proposed as a preliminary 
prototype for system evaluation. 

A logic design was formulated, and its veracity was 
analysed via a simulation using actual sensor word lengths 
and appropriate fixed point arithmetic. The simulation as 
a design tool verified an optimum design with proper 
trade-offs between system performance and power and size 


projections. 
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II. DESIGN CONSIDERATIONS 


A. SYSTEM EQUATIONS 
Although the four input variables of temperature, 
pressure, and the two current components are sampled, only 
the pressure data is sampled for unusual events in the 
tsunami spectrum of 1 to 30 CPH [Ref. 1]. A more sophisti- 
cated device could use any number of detectors on any of the 
inputs using the same temporary storage buffer memory. 
1. First Difference and Limiter 
In order to remove the mean pressure and long-term 
drift, the first stage in the detector is the first dif- 
ferencer. It computes the difference between successive 
oles Peo - Pe]. in Order to pretect the detector -rem 
any spurious data which would give a false indication, a 
limiter which limits the first difference according to the 
following relation is the next step in the detector: 
Eel 2a 1 Eee 
Ne = 


Sg esa) Lt Oe Sas (61 


Pee recor i iter 
The tide filter is a fourth order filter with zeroes 
chosen to remove the diurnal and semidiurnal tides. Two 


Sections of the form 


Sy a5 a1 Sea + bi 54-2, = feceor) Wy (r) ae where 


W,(t) = 1 + a,6(d) + by, 6 (24) 


ae 





with aj, by and a>, bog chosen to give the required zeroes 
are combined by convolution to give the fourth order 
difference equation: 


W(z) = 1 + azaz6(A) +[ayag + byb9]6 (24) + 


[azb2 + a2b1/6 (34) + by b56 (44) 
With A = 1/2 hour, the five inputs to the filter will be 
the first difference SI chs the times 0, -30, -60, -90, and 


Seeweminutes,. The constants are calculated to be 


ag = iL 
eee ag = aoc wee 
aya ate bi b> = 5.841800 
a,b, cob = 7S. aeBo2 

eYaey 

oo 6 inkegiration 


To obtain the power spectral density of the tide 
filter the signal is next Squared and averaged over a fading 


Function with a 60 minute time constant. 
P2(t) = aPZ(t-A) + P2(t) az .991666 


The only function remaining is to determine if the power has 
exceeded a preset threshold indicating significant energy in 


the tsunami band to signal an "event." 


a Contecol FuncElons 
The system proposal calls for a dump of the previous 


two hours of one-minute samples whenever a tSunami is 
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detected. This will ensure that an adequate history 

iooing Up .tO.an event is stored. Alsomineluded in the 

Hic is the circuitry to continue storing the one-minute 
samples until the second consecutive two-hour interval which 
is below threshold. Extra logic gates are provided, however, 


pemecevariety of controls lending flexibility to the 


detector as a "black-box" component. 


B. SYSTEM ARCHITECTURE 

The logic design for a special purpose digital machine 
begins with the system equations, block diagrams and a gen- 
eral understanding of overall system requirements. Special 
purpose structures can usefully be described by the general 
system pictured in Figure 2 illustrating the two major parts 
of a special purpose machine, the data processing and 
timing and mode circuitry. 


So lny 
Ne US 













DIGITAL 
PROCESSING 
CLREURER TG 


CIRCUITRY -————— 





SYSTEM 
OUTL UTS 


Figure 2. Special Purpose Structure 


13 





The data processing circuitry includes all input and 
output interfacing circuits, the data operators and system 
blocks, any storage of memory devices, and, in general, the 
necessary hardware to take system inputs and derive system 
outputs. The timing and mode circuitry serves much the 
Same purpose for the special purpose structure as the soft- 
ware program does for the general purpose computer. The 
timing circuitry sequences operations in a set pattern from 
some initial time. The mode circuitry controls operations 
which are squalene of the data being handled as well as 
time. 

In considering the best approach to the design of the 
digital package, the need for low power has repeatedly been 
stressed. For this reason, serial processing circuitry 
which has the advantage of economy of hardware at the ex- 
pense of speed is used predominantly with only slight 
modifications to meet other of the specified system 
requirements. 

1. Memory 

The major component and power requirement of the 
detector will be the temporary data storage section or buf- 
fer memory. Investigation into this problem centered around 
core memories, shift registers, thin films and delay lines. 
Core memories have the definite advantages of nonvolatility, 
simplicity, zero quiescent power, reliability and random 
access capability. The same is basically true of the new 


thin film adaptations of* ferrite storage. Shift registers 
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are popular as short term buffer storage, are serial in 
Operation, and require no external read ee Reavare. 

Since random access, nonvolatility of information, 
and access time are of little importance here, the most 
favorable devices are long shift registers which are compa- 
tible with the logic family. A discrete state shift register 
of length 10,000 bits at first seems to be rather prohibitive, 
but recent trends in MSI are making these long shift regis- 
ters well-suited for small size memories. Using current low 
power MOS technology as many as 64 static bits are available 
on a Single chip. Since such long shift registers used as a 
recirculation memory in a serial machine need Simply a 
Single port input, a single port output and a clock, the 
total amount of hardware involved for the complete memory is 
very small. 

Also considered were dynamic shift registers, which 
have more than three times the number of bits per chip as 
BiemerG Gbegistera (Ref. 2]. Although the quiescent power is 
lower, the total power is greater than for static devices 
because of the transient power of continuous clocking. The 
resulting long static discrete state shift register chosen, 
therefore, was dictated by available hardware at low power, 
serial operation for simplicity, recirculation at an extreme- 
ly low rate for reduced hardware count and data accessibility, 
and low clock rate for low duty-cycle. The organization of 


data in the shift register memory is shown in Figure 3. 


ps 





2. Data Organization 


The type of data organization and corresponding 
arithmetic code to be used for the detector results from 
the form of the input data. All five input channels (pres- 
sure, temperature, current components and time) will be 
Sixteen bit binary magnitudes. For this reason, the Data 
can be simply word organized within the samples which are 
sequentially organized as shown in Figure 3. The arithmetic 
code chosen was initially a one's complement code for simpli- 
city. Further study, however, indicated the use of a two's 
complement code waS more compatible with serial processing 


and would ultimately require fewer operations and hardware. 








S6SO BITS 





TIME CURRENT, CURRENT,, TEMPERATURE PRESSURE 


Figure 3. Data Organization in Memory 
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Two additional considerations which will be clarified in 
the later section regarding the use of RCA's COS/MOS logic 
family are the use of positive AND/OR logic for initial logic 
design and the use of asychronus timing circuits. The short- 
cuts involving the use and marriage of NAND and NOR logic to 
reduce the package count is a matter of circuit juggling to 
make optimum use of the available gating packages. Asychro- 
nus or pseudosynchronus timing of all operations again 
results from the tradeoff of speed for simplicity. 

Be. Flexibility F 

Two of the specifications for the unusual event 
detector are compatibility and flexibility. The input sec- 
tion of the device is therefore designed to accept a variety 
of input combinations. ‘The digital sensors for the deep sea 
in most cases are frequency modulated signals corresponding 
to the particular input variable. They are measured by 
allowing the output frequency to count a binary counter for 
some predetermined period [Ref. 3]. 

The input section is designed to accept five 
variable frequency or variable period signals. The exact 
period of measurement (from one second to four minutes) or 
the choice between frequency or time measurement is easily 
changed, using removable solder jumpers to program the input 
in a variety of fashions. 

The storage of the 120 minute samples in a static 


register provides a flexibility of output. A simple 
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clocking gate or even an external clock can be used to read 


out this information in almost any format. 


C. IMPLEMENTATION WITH COS/MOS 

Recent trends in semiconductor technology toward medium 
and large scale integration have made the MOS transistor 
increasingly attractive due to its ease of manufacture and 
versatility of operation. Single chip building blocks of 
rather large arrays, registers, and memories are being incor- 
porated into many logic families; of primary interest here, 
however, is the use of the MOS device in a complementary 
symmetry arrangement which has the distinct advantage of 
extremely low power [Ref. 4]. 

At present the leading supplier of a complete logic line 
using complementary symmetry MOS transistors 1S Radio Corpo- 
ration of America with its COS/MOS line. The COS/MOS CD4000A 
series of integrated circuits features a variety of options 
including special low voltage types, high reliability cir- 
cuits, various packaging styles, and the usual military or 
industrial temperature ranges. 

Pe) Cone lenientary so ynneeny (el aed tiny 

The breakthrough that made complementary symmetry 
logic practical in medium scale integration was the ability 
to fabricate opposite polarity, but Similar devices in large 
numbers on the same chip. The basic gate circuit from which 
all other logic functions are derived is shown in Figure 4. 
The usual load resistor is replaced with another FET of 


opposite polarity; thus, there is no quiescent current 
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except for the leakage currents on the order of nanoamps. 
The pair of MOS transistors are simply switches which tie 
the output to either the supply voltage or to ground. In 
true complementary logic, the circuit for the. P devices is 


the dual network function of the N units. 





fe) 
Vo 
INVERTER 
N 
ite N 
N 


® IN d IN b IN 165 


Figure 4, Three Input Complementary Symmetry MOS NOR GATE 


The advantages derived by the complementary 


arrangement are more than just the microwatt standby power: 


1a, 





a. Simple one supply direct-coupled logic using 
only MOS transistors 

b. Large fanout capability because of capacitive 
inputs 

c. Good noise immunity due to large threshold 
voltages 

ad. High speed operation because output node 
Capacitance is always charging or discharging 


through the "on" unit 


The P and N channel units can form a shift register 
using true complementary logic with two flip-flops per bit 
for two phase operation; however, a Significant reduction in 
components can be obtained by the use of (1) the excellent 
Capacitance storage feature of the MOS gates, and (2) the 
two halves of the MOS flip-flop, one of which will provide 
the function of delay [Ref. 4]. As shown in Figure 5, 
Opening the switches S, and S5 which form the Cross -—COoup@ing 
connections results in two separate storage elements, each 
capable of storing for an extended but finite period of 
time. In this state, new information can be input to the 
left gate, while the level of the right half is transferred 
to the next stage on the right. Closing the switch S,5 and 
then closing S, locks the one-zero information indefinitely. 
The two phase clocking requires a rather quick transfer of 
data from the left to the right halves of a stage, but an 


indefinite delay is allowed between stages. 
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DATA 
TRANSFER 


Figure 5. Basic Shift Register Stage 


3. Power Dissipation 

Considering the total power dissipation as the sum 
of both the quiescent and transient powers, further power 
reduction can be obtained by proper design of the switching 
characteristics. Minimizing the logic level swings, using 
true complementary loads, and minimizing leakage currents as 
well as all output and load capacitances form the general 
approach to an optimum system. The total power P dissipated 
in a complementary set-reset flip-flop is given by 

Pag’ +P. 
where Co is the node capacitance, Vo is the supply voltage, ~ 


f is the operating frequency, and P, is the standby power. 
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For low switching power, the supply voltage Vg must be 
poaweed to the lowest possible value. AC a Vg Of -4 volts, 
and an output capacitance Co of 23 pf the frequency depen- 
dence is shown in Figure 6. Note, also, the large increase ~ 


in quiescent power with temperature [Ref. 4]. 


MICROWATTS 
FOR A 

SINGLE 10 
Hore =-FLOP 
STAGE 





) s) 10 IRS 20) KZ 
Figure 6. Power Dissipation as a Function of Repetition Rate 


Since only a fixed number of operations are completed 
each minute in the detector, the clock rate will have no ef- 
fect. The power economy is gained by reducing the total num- 
ber of redundant operations. The repetition rate can go 
down to DC; however, the longest allowable high state of the 
Clock input of the 64 bit shift registers can be no more than 
100 microseconds. For this reason, the clock frequency was 
chosen to be 19.360 kHz which provides a good margin of 


safety. 
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III. HARDWARE IMPLEMENTATION 


The logic design formulated in this section is organized 
according to the logic blocks diagrammed in Figure 1 of the 
Introduction. All logic functions are described and defined 
in the text and the figures. The later sections appraise 
the cost, Size and power requirements of the resulting logic 


packages. 


A. CLOCK AND TIMING 

The first section of the detector, which sequences all 
operations and controls the functions of the device, is the 
clock and timing circuitry. Table I, which may be consider- 
ed as the program, lists all sequences and defines the 
timing signals originating in the COS/MOS counter circuitry. 

Since the COS/MOS inverters have nanowatt standby power 
and microwatt operational power, an excellent low power clock 
can be designed using the forward gain of an inverter with a 
Suitable feedback arrangement. A good amplifier section for 
the oscillator is formed with an inverter stage and feedback 
resistor R, large enough (ePOMmeorthal the eaeremma one 
phase of the feedback are not appreciably affected. The 


resistor biases the output at about half the supply voltage. 
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R 22 k 311.6 kHz 


220 pf 8-50 pf [47 of 


Pigune 7. Crystal @seillacor Cireulte Using 
RCA CD4007 Inverter 

The crystal controlled phase shifting network is added 
in Figure 7 with provision for trimming to the design fre- 
quency. Using this circuit, Stability should be on the 
order of 5 ppm and total supply current from a supply voltage 
of 4 volts is 50 microamps [Ref. 6]. 

Frequency division of the crystal oscillator is 
necessary due to duty cycle considerations which place the 
logic clocking frequency at 19,360 Hz as shown in Figure 8. 
One CD4004 provides four flip-flops in an asychronus counter 
to divide the 311.6 kHz crystal frequency down to the clock- 
ing frequency of 19,360 Hz. The bit counter consists of one- 
half of a CD4013 dual flip-flop plus a CD4022 eight bit 
divider/decoder to aid in decoding the timing signals. The 
word counter is ones CD4022 connected to divide by five 


with a CD4001 gating package to reset the counter. 
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The sample counter is a CD4004 seven bit counter 
connected to count to 121, and ten other logic packages 
provide the necessary reset and decode functions. The 
second counter is a similar arrangement connected to reset 
at 60 using the same type logic packages. The first stage 
in the second counter is the recirculation counter; it 


divides by two to give a frequency of 1 Hz. 
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TABLE I. TIMING 


mnese Boolean expressions define sequentially the timing 
signals referred to on the later diagrams. The Gating “Cimes 
cuitry for these expressions is shown in Figures 10,1l1l,and 12. 


Sec(0)-R(0) RECIRCULATE BUFFER MEMORY 
Sec (0)e¢R(0) FIRST DIFFERENCE AND LIMITER 
S(L)*Ep = Sec(0)*R(0) [S(0) +S (30) +S (60) +8 (90) +5 (120) | 
S(H)*E, = Sec(0)*R(0) [S(1)+S (31) +S (61) +S (91) +S (121) | 

Ep? W (0) Clock 16 bit delay register 
and adder 
S(H)«Ep[W(0) + w(2)] Shift Limiter register 
S (H) e Ep? W (1) Check difference and Set 
Takecn, LE gGocd.) "Sees ten 
change latch if necessary 
S (H)e Eps W(2) ; Shift present difference 
EO EGNPOLranyY Steemade. 
S(H}eH,s W(2)°B(15) Change sign. 
Sec(1)+R(0) TEMPORARY STORAGE 
EpeW(2) [B(8) to B(15)]* Shift in 8 bit difference 
Esum’ W (0) Set t OUmmett ference 
Sec (2)*R(0) DEE Lie Re 


Esum = Sec(2)¢R(0) 
Esum |W (0) B(5) toB (15)]*+W(1) Multiply by -3.920392 
E sum |W (0) B(7) toB(15)}*+W(1) Multiply by 5.841800 


Esun' S$ (1) + $(121)] » Add differences 1 and 121 
(1)*B(4) to B(12)]* 


E sunt [S (1) +S (121)]}ow(1) +B (12) Set sign latch 


1B sea + [Esq * W(1)+B (0) Shaert accumulator 
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TABLE I. Continued 


Sec(3) R(0) SQUARE 
Egq = Sec(3)+R(0) [S(0) to $(75)]* 
Esg*W(1)* B(0) Shift accumulator in’ sum 
and checke® bit for "one" 
Esq? W (0) Shift input store and SQUARE 
accumulator 
Sec(4)e«R(0) INTEGRATE 
Brn ame’ (4) 5 Ri 0}e S(O) 


Eint W(0) [B(11) to B(15)| Shift in the squared value 


Bee) Shift registers 


Sec (59) *R(0) INPUT SECTION 
Ein = Sec(59)¢R(0)*¢S (121) 
EsneW (0) Shift in pressure value 


Eine W (1) Shift in temperature, 
reset pressure counter 


Eine W(2) SJaah ime, Ebot ceibieieeone 57 
reset temperature 


Ene Wis) Sheet. a current, 
reset current, 

Ein? W(4) Shift in time, reset 
current,, 

7B (a) Seo Big) =| Bi) + SBiatl +... + Bia) 


The decoding circuitry for these signals are shown 
together in Figures 10,11, and 12. Usually they will be 


located on the samé board as the device they control. 
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Bee EROCESSING CIRCULTS 
L, @enpaut Section 

The five inputs consisting of variable frequency 
Signals are used to either count or gate five sixteen bit 
binary counters. The fifth counter is provided with the 
eaeidativonal circuitry so that it will count time in minutes. 
On command from the system clock each of the five counters 
sequentially shifts its value into the system memory and, 
S<ecpt for the time channel, resets £or the next counting 
interval. Figure 13 shows only the pressure and time channel 
because all five inputs are the same except for the gating 
on the minute counter. Since the two-hour sequencing signal 
for the detector is obtained from the time counter, it must 
reset at 65,520, a multiple of 120 minutes, so that the 
two-hour sequence iS maintained. 

All of the counters are made from CD4004's. Three 
counters for each input provide the necesSary sixteen bits. 
Two eight bit parallel to serial shift registers convert 
the counter value into serial form for the system's buffer 
memory. Fourteen additional packages provide the gating and 


mMO@e  COntCrOl Circuitry . 
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25 System Memory 


The system memory consists of a discrete state static 
shift register of 9680 bits in a recirculation loop storing 
the past 121 minute samples organized sequentially by words 
and bits. The bit, word and sample counters provide for 
addressing the data; further division of the clock frequency 
into seconds and minutes provides all necessary timing sig- 
nals for the rest of the detector. The timing signals are ob- 
tained by coincidence gating across the stages of the asychro- 
nus binary counters. In order to reduce the transient power, 
the data is recirculated through the register length only once 
per minute. On this one pass, clocking and gating read out 
the ten pressure words needed for computation of the five 
first differences. The 9680 shifts averaged over a minute's 
time results in an average clocking frequency of 150 Hz. 

From data Similar to the graph of Figure 6, the transient 
power used in the memory should be about 100 microwatts; 
however, the 64 bit shift registers are still in the devel- 
opmental stage and this power estimate 1S meant only as a 
comparative estimate. 

In order to obtain the 9680 bit length in COS/MOS 
logic, 152 logic chips must be series connected. The exact 
9680 bits required is made of 151 sixty-four bit TA5989 
registers and one CD4006 sixteen bit register. The total 
memory package should dissipate about 0.3 mW and fill a 
volume of 0.34 cUu.tteerwon six Stamaard 2.5 by 5.0 ianech Logie 


cards. 
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3. First Difference and Limiter 

First differences are calculated once a minute on 
the five sample pairs S(0), S(1); S(30), S(31);...S(120), 
S(121) by delaying the first sample in a register, comple- 
menting the next and adding the two serially. Since the ori- 
ginal input number is aintenes to overflow the input counters, 
the signed result must be corrected if the sign is inconsis- 
tent. Gating also checks the difference for unreasonable 
values before allowing them to pass on to temporary storage. 

The sixteen bit delay and all the eight bit storage 
requirements in the first differencer can be Beanies fe Tae Th) 
various connections of the CD4006, an eighteen stage regis- 
ter. Two CD4015A dual four stage parallel output registers 
provide the necessary connections used in the limiter to 
sense mrcorrect the difference value. One-third of a 
TA5963 triple serial adder is used for the adder anda 
CD4013 dual flip-flop performs the complement function. 
Twelve additional packages as listed in Table II provide 
all the gating and control logic. 

4. Temporary Storage 

Since only correct first differences are allowed 
through the limiter, the temporary storage must be capable of 
nondestructive read-out when the tide filter calls for inputs. 
This 1s a very simple matter of circulariy shifting the data 
to the right. Only the 8 most insignificant bits need to be 
stored, because the average first difference is on the order 
of 4 bits. The validity of this assumption was verified by 


Simulation. 
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By tice Filter 

The five differences temporarily stored are simply 
Multiplied by appropriate weighting values and summed to ob- 
tain the filtered output. Since there are only three dif- 
ferent constants, the differences are gated to the proper 
section through serial switches and the sum of these values 
is stored in a register which also serves as the input 
register to the squarer. 

Since one of the weighting constants is unity, only 
two multipliers are required. They simply shift and add the 
input the appropriate number of times and proper clocking 
keeps the result in the same sixteen bit format. 

The constant -3.920392 can be approximated with the 
number -3.920410 or 22 Fo a + 27-6 42-10 + 2-11). This 
constant, when multiplied by the input number of sixteen bits 
gives a possible maximum length result of 29 bits. Through 
Simulation and device consideration, however, the input 
difference will rarely exceed eight bits and the approxima- 
tion to reality of the constants makes the significance of 
the lower five bits questionable. Truncation to sixteen 
bits before addition, then, will have no effect. Similarly, 
5.84180 is approximated by 5.842797 or 23 ~ (21 + 2-3 + 275 
+ 2-9). 

The maximum of the sum of the five weighted differ- 
ences will be approximately the sum of the magnitudes (1 + 4 
+ 6 + 4 + 1) times the average input value. Limiting the 


maximum first difference to seven bits or 128), requires 
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eleven whole number bits in the accumulator leaving four 
Beactional bits plus the sign bit. To multiply by the con- 
Stant one, the input need only be left-shifted four bits to 
account for the change in the binary point. Clocking and 
gating accomplish the necessary four bit delay. 

The analytic two's complement code used throughout 
the detector requires no special handling in the multipliers; 
however, Since the eight bit differences are to be extended 
to sixteen bit precision words in the accumulator, care must 
be taken to pack the additional eight created bits with 
zeros or ones depending on the sign of the original number. 
The sign latch in Figure 18 tests the sign from the lst and 
l2lst difference and ensures that the input is held at the 
proper level if there are no incoming bits. 

All of the timing signals referenced on the figures 
are defined in Table I. The Figures 10, 11 and 12 diagram 
how the timing signals are derived through gating and 
decoding the frequency divider outputs. 

The total number of twenty-five delays in the two 
multipliers is provided by four CD4015 dual four bit regis- 
ters. Three more TA5963 triple serial adders perform the 
summing functions, with the result stored in an additional 
two CD4015's. The two complementors are similar to the one 


used in the first difference and diagrammed in Appendix A. 
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6. Squaring 

The squaring package observes the value stored in 
the sum accumulator of the filter and shifts and adds it into 
the output storage. The square circuitry would handle nega- 
tive as well as positive numbers if care were taken to pro- 
vide all the necessary ones in the extra negative bits. A 
Simpler method to implement, however, is to simply rectify 
the input. By testing the sign bit in the sum accumulator 
before each operation, negative numbers can be complemented 
bit by bit. This method, as shown in Figure 19,,gives an 
error of one in the least significant bit because it takes 
the one's instead of two's complement. The error is 
insignificant. 

For each bit the input is shifted and if the current 
bit 1s al, the shifted value is added to the result. The 
sum in the output of the filter contains a sign and 0) Dg 
fractional bits; in order to make the detector adaptable to a 
variety of dynamic ranges, the tap on this register is made 
adjustable. The input to the squarer can only be eight bits 
in length, but the adjustable tap gives a range of the most 
sensitive eight bits from 1/32 to 15 to the maximum values 
between 16 and 65,536. Simulation later showed the most 


probable tap to be at the binary point. 
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7. Integrator 

Since the squared output is to be integrated over a 
decaying time constant of sixty minutes, the maximum avail- 
able storage must be about thirty times the maximum value to 
ensure against overflow. To stay within the sixteen bit pre- 
cision, the input squared value must be truncated to the 
eleven most significant bits. The decay factor of .991666 
is accomplished by delaying seven bits and then subtracting. 
Taps are accessible on the delay register to vary the time 
constant of integration. Coincidence gating on.the storage 
register synchronized to the clock determines the value in 
the tsunami detector and provides the control logic with its 
go/no-go Signal. 

The necessary storage for the squarer and integrator 
can all be formed from the CD4006 eighteen stage register. 
A CD4015 dual four bit register provides for the Minted > Ga — 
tion in the integrator along with two of the summers left 
over from the triple serial adder used in the first differ- 
encer. The complementor again comprises a dual flip-flop 
CD4013 and nine additional gating packages complete the 


Squarer and integrator. 
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8. storage Control 


This section of the detector is designed with flexi- 
bility as the prime consideration since the actual device used 
for long term storage will be an independent system. Whether 
a tape system or core storage is used the only output required 
of the detector is a yes or no answer every two hours quali- 
fied by considerations of its previous answers. A critical,or 
event, is defined whenever the value in the Integrator regis- 
ter exceeds the threshold set by the sense tap. A store 
command is given under the following conditions: 

a. A critical is detected and the past n consecutive 

hours have not been critical; 

b. A critical is not detected and the past two hours 


were critical. 


UNUSUAL 
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RO RO : 
> N BIT COUNTER 
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Sec(0) MIN (120) 


Figure 21. Control - Three Flip-Flops and a Counter 
Provide Logic to Determine if an Unusual 
Event Should Be Signaled 
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C. PACKAGE CONSIDERATIONS 

The COS/MOS line offers sufficient complexity per chip 
to keep the overall package count low, and also, enough var- 
iety to provide for optimum utilization of the available 
circuitry. The total package requirements for each of the 
sections of the detector are listed in Table II, with a 
cost estimation in Table III. 

1. Space and Cost 

The size and cost predictions of the detector are 
based on an experimental set-up which will use dual in-line 
packaging on a general purpose bread board arrangement. 
Standard Logic Inc. offers a line of logic cards which have 
various combinations of fourteen and sixteen pin DIP sockets 
and all mounting and interconnecting hardware. Computer aided 
Wiring and interconnections complete the system. Using the 
Standard Logic boards from 20 to 30 DIP packages depending 
on the number of sixteen and fourteen pin packages, can be 
arranged on each 2.5 by 5.0 inch card. Figure 22 shows the 
logic arrangement and overall dimensions of this experimental 
set-up. 

For reliability the detector should be made with 
flat-packs on printed circuit boards with all solder con- 
nections. The space requirements will be approximately 40 
to 60 percent of the experimental arrangement due to the 
elimination of the DIP sockets and the interconnecting 


wiring. 
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TABLE II. PACKAGE REQUIREMENTS 


FUNCTION FIGURE NUMBER LOGIC CIRCUIT 


CLOCK 8 2 


CD4022 
CD4004 
CD4013 
CD4009 
CD4023 
CD4001 
CD4009 


CD4004 


CD4013 


INPUT ILS) ibs) 


10 


eee 


ee, 
MEMORY 14 sBSy 


2 


ts 


154 


fe @erciel. 
CD4004 
CD4014 
CD4011 
TASS87 2 
CD4023 
CD4000 
CD4013 
CD4009 
BOrtat 
TA59 89 
CD4011 
CD4006 


total 


50 


Divide by 8 Counter/Decoder 
Seven Stage Counter 

Dua Fl 1ip-f£ftep 

Hex Inverter 

Triple 3-Input NAND Gate 
Quad 2-Input NOR Gate 

Hex Inverters 

Seven Stage Ripple Counters 


Die Fre fee 


Seven Stage Ripple Counter 
8-stage Parallel in/Serial out 
Quad 2Z2-Input NAND Gate 

Dual elite -f£lop 

Triple 3-Ghput NAS) Gaze 

Dual S-input NOR ois tne aecns 
DiS IL 1D) Seba h sy—18 Ihre)», 


Hex Inverter 


64~Stage Static Shift Register 
Quad 2-Input NAND Gate 


18 stage Static Shift Register 





TASLE IL. Continued 
FUNCTION FIGURE NUMBER 


FIRST 15 5 
DIFFERENCE 


ao 


LOGEC CikGgeT 


CD4011 
CD4013 
CD4015 
CD4016 
CD4012 
CD4009 


CD4006 


TA5963 


17 


TEMPORARY 16 3 
STORE 


total 


CD4006 
CD4011 
CD4016 


CD4009 


CD4023 


MULTIPLIERS 17 4 


total 


CD4015 
TA5963 
CD4011 


CD4013 


CD4023 


13 


total 


5) 


Quad 2-Input NAND Gate 

Dual D Flip-flop 

Dual 4-Stage Parallel out 

Quad Bilateral Switch 

Dual 4-Input NAND eae 

Hex inverter 

18 Stage Static Shift Register 


Triple Serial Adder 


18-Stage Static Register 
Quad 2-Input NAND Gate 
Quad Bilateral Switch 
Hex Inverter 


Triple 3-Input NAND Gate 


Dual 4-Stage Parallel out 
Triple Serial Adders 

Quad 2-Input NAND Gate 
Dual D Flip-flop 


Triple 3-Input NAND Gate 





TABLE II. Continued 
FUNCTION FIGURE NUMBER 


SUM 18 2 
ACCUMULATOR 


f= 


Fe 


BHOGIC CRPRCUIT 


CD4011 
CD4015 
TA5963 
CD4023 


CD4000 


CD4002 


© 


SQUARE IS, 2 


SO ie 


CD4006 
CD4011 
CD4001 
CD4009 


CD4016 


TASI6S 


INTEGRATE 20 3 


ib 


_—L__ 


i 


total 


CD4015 
CD4011 
CD4013 
TA5963 


total 


S72 


Quad 2-Input NAND Gate 

Dual 4-Stage Parallel out 
Triple Serial Adder 

Triple 3-Input NAND Gate 

Dual 3-Input NOR plus Inverter 


Dual 4-Input NOR Gate 


18-Stage Static Register 
Quad 2-Input NAND Gate 
Quad 2-Input NOR Gate 
Hex Inverter 

Quac Bilateral Switch 


Triple Serial Adder 


Dual 4-Stage Parallel out 
Quad 2-Input NAND Gate 
Duals PP vip. Lep 


Triple Serial Adder 





TABLE II. Continued 


FUNCTION FIGURE NUMBER 


EOGiC CLEReury 





STORAGE 2 il 2 cpd4011 

CONTROL 
1 cpd4004 
LPep2013 
1 +‘@b4009 
5 total 

TIMING 10 6 “CDdaot1 

GATES ‘lg 

& 12 2 ~CD4001 

2 eb40 24 
SD OS 
1 cpDd4002 
1 ecp40l2 
2 CpD4009 
16 total 


23 


Quad 2-Input NAND Gate 
Seven Stage Ripple Counter 
Dual Flip-flop 


Hex Inverter 


Quad 2-Input NAND Gate 
Quad 2-Input NOR Gate 
Triple 3-Input NAND Gate 
Triple 3-Input NOR Gate 
Dual 4-Input NOR Gate 
Dual 4-Input NAND Gate 


Hex Inverter 





STWNOIS ONIWIL 
LOAd LO 

TOY LNOO 

SLOANI 

YaMOd 


YOLIANNOO 


TOU LNOO 
A LVYOudLNI 


uoTzeanhHtyuoy TeotsAud 


TOULNOO 
GQNY FLVYOULNI 
aav~NOs 


qaWo T SdId Sd 


YOLVTAWNOOV 
WOS ANY 
odd TIdI  LTAWN 


quwo T 


pala 





°77Z oanbtd 


qOVY OLS 
AUYVHYOdWNdL UNV 
GONGYETATd LST 


quvo T paid Os 


Sian 9) oe oda 7€ SL AY OWA 
Sadao c sdid Sv LOdNI 
quvo T SdId O€ MOOTO 


YOLWINWNDOW Pry \ gyoOLs aNw = 7 = 


CNW 


Sad TIdi LiInw AONGYAATATG Loar 


AY OW YIN LOdNT MOOTO 


54 





TABLE III. COST EVALUATION 


Low ‘Yoltage Series CD4000A 3 to 15 volts 


Flat Pack =5or Geto. Zac? 
Number Cost Each. Total Cost 

CD4000AK 2 $ 7.00 $ 14.00 
‘CD4001AK 6 7 AN 44.40 
CD4004AK 19 ee, 7s 280.25 
CD4006AK 7 18.00 26 G0 
CD4009AK 12 11.00 132.00 
CD4011AK eu me. 40 22.260 
CD4012AK 3 7265 22.95 
CD4013AK 9 10) .0)!5 90.45 
CD4014AK 10 Ao) Ls AOL 0 
CD4015AK 1d 720) Pa 222275 
CD4016AK 6 10.80 64.80 
CD4022AK 2 20.95 41.80 
CD4023AK 19 7.65 Las: ails 
CD4025AK 2 7.65 15) 80) 
A587 2 ioc ove 30.00 
TA5963 7 14.00 100.00 
TA5989 151 20.00 3020.00 
ROMAT, $4752.55 


Similar arrangement in 


Dual In Line Plastic Packages Approx- $2000.00 


—ACreror + 15 C 


SS 





2. Power 

The actual quiescent power required varies 
considerably from chip to chip and is very temperature 
dependent. In most cases, RCA specifications include a 
"typical" and "maximum" value. The estimates given assume 
the "typical" measurement and an ambient temperature of 
25°C. The actual power at the ocean floor, therefore, 
should be less than the estimate. Also the transient power 
is of no significance in the digital processing circuitry 
because the maximum number of operations of any of the cir- 
cuits is still on the order of 160 per minute. This gives 
an average repetition rate of only 3 Hz when averaged over 
the one minute period. For this reason only the transient 
power of the clock, the decoder and the buffer memory are 
included. 

The total power drain of 0.5 mW required for the 
whole detector system is almost three orders of magnitude 
smaller than a comparable logic package of conventional cir- 
cuitry. Compared to the power requirements of the sensors, 
accoustic links, and other functions of the capsule, the 
power for the detector is negligible. The power for one 
year's operation should be obtainable from a single 2 amp- 
hour storage battery. The power used to run the crystal 
oscillator should be well regulated and be run from a sepa- 
rate mercury battery for voltage stability. All of the other 
detector circuits can be supplied from lead-acid cells due to 


the excellent voltage compatibility of the COS/MOS line. 
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If power drain were not so important, more than a 
two to one reduction in package count could be effected by 
rebuilding the memory out of conventional MOS registers, 
The power required would increase almost ten-fold. Future 
developments in the COS/MOS line toward LSI might even 


reduce the size as much with no sacrifice in power. 
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TABLE IV. 


LOGIC CURRENT 


CD4000 
CD4001 
CD4002 
CD4004 
CD4006 
CD4009 
CD4011 
CD4012 
CD4013 
CD4014 
CD4015 
CD4016 
CD4022 
CD4023 
CD4025 
TA5872 
TA5963 


TAS9 89 


Dual 3-Input plus Inverter 
Quad 4-Input NOR Gate 

Dual 4-Input NOR Gate 

Seven Stage Ripple Counter 

18 Stage Static Register 

Hex Inverter 

Quad 2-Input NAND Gate 

Dual 4-Input NAND Gate 

Diabet) wh hohe 

8 Bit Parrallel in/Serial out 
Dual 4-Stage Parallel out 
Quad Bilateral Switch 

Divide by 8 Counter/Decoder 
Triple 3-Input NAND Gate 
Triple 3-Input NOR Gate 
bocive gd) K bli p= hep 
Triple Serial Adder 


64-Bit Static Shift Register 


TOTAL STATIC POWER 


TRANSIENT POWER: 


CLOCK 


MEMORY 


TOTAL TRANSIENT POWER 


aye 


POWER REQUIREMENTS 


STATIC POWER NUMBER 
Microwatts 
Om Onl: 2 
OO 6 
oF 2 
i ibs, 
OR05 7 
C205 12 
0.005 Sik 
C5 OOS, 3 
Om025 2 
4.0 10 
4.0 il 
OAL | 6 
lgels 2 
0.005 19 
OmOe 2 
On Oz 4 
eo 7 
0 ILL 


268.8 Microwatts 
270 
110 


380 Microwatts 





IV. COMPUTER SIMULATION 


A preliminary step in the design of the system components 
was the evaluation of the required functions through computer 
Simulation. The general algorithm has previously been ana- 
lyzed [Refs. 1,5]. System evaluation of the actual logic sys- 
tem should be as close to actual hardware circuitry as 
possible with fixed point numbers and identical computer logic 


Simulation. 


eee 6COMCOR CI 5000 

For actual checkout of some of the functions in hardware, 
mewatgital patching section of the CI 5000 which us located 
in the Electrical Engineering Computer Lab, was used. The 
complementor for the two's complement arithmetic was the 
first function wired. The second logic block that was test- 
ed extensively in hardware before design was finalized was 
the multiplying scheme. Since a single multiplier used all 
32 of the computer's general purpose flip-flops and a major 
portion of the NAND gates, this method was applicable to only 
the system blocks; another computer had to be chosen for the 


overall system evaluation. 


me CDC Lod 
With the consideration that the simulation should retain 
as much simularity to the actual hardware as possible, the 


CDC 160 computer, also located in the Electrical 
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Engineering Computer Lab, was chosen, The distinct 
disadvantage, later to become obvious, was the twelve bit 
word limitation in trying to check out a sixteen bit system. 
In most areas of the program and logic flow, there were no 
differences between the designed sixteen bit system and the 
twelve bit simulation; however, the checkout of the multi- 
pliers required further investigation. At this point the 
168 arithmatic unit was used as peripheral equipment to 
allow for double precision. 

i” Adaptation to Frxed Point 

A magnetic tape wasS available from the 1.G.P.P. 

containing two records of pressure variations at the 2000 
fathom level in an area near San Clemente. The data had 
been taken at four minute intervals over a four day period 
with both a Vibratron and a Hewlet-Packard sensor. The data 
was linearly packed to simulate the required one minute sam- 
pling rate, and this prepared tape was used [Ref. 5]. A 
major portion of the program listed in the Appendix was re- 
quired to read the BCD tape and convert it into pure binary 
form. Several problems due to word length limitations were 
solved by reading in only the four least significant decimal 
digits and converting them to binary through a table look-up 
scheme. Simple testing of the resultant value for consis- 
tency removed any ambiguity in operations with the numbers 
except when a change in the fifth BCD digit occurred as noted 
on Figures 24 and 27. This method was incorporated into the 


actual design also because the input counters were allowed 
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to overflow any number of times. Figures 23 and 26 show a 
record of the tape data compared to the truncated version 
used on the 160. 

2. Precision Considerations 

In order to avoid going to double precision and using 
the external arithmetic device, which would no longer bear 
any, relation to the hardware being simulated, careful study 
of the data revealed several compromises. In most cases a few 
of the most significant bits could be dropped because the 
particular data did not reach the maximum excursion designed 
for in the sixteen bit system. This, combined with truncat- 
ing one or two of the less significant bits, reduced the word 
length to the required twelve bit format. No significant 
degradation in overall performance was noted from these com- 
promises except in the case of the multipliers. Here, twelve 
bit precision allows for an error of .002 which is not insig- 
nificant since the least-count resolution of the input sen- 
sors is on the order of .0005. The multipliers were tested 
mole proper functioning in twelve point individually, but, te 
allow for system evaluation and to obtain usable results, 
double precision was used. 

The only other notable difference between the 
Simulation and the design was the use of one's as opposed to 
two's complement arithmetic. Although one's complement re- 
quires the additional operation of end-around carries to 
correct the sum on subtraction, this is handled in the ma- 


chine; therefore, on the software level, compensations were 
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made so that the simulation was the same as the two's 
complement on the hardware level. 
3. Simulation Results 

The program operated on the input data samples in the 
same format as the proposed system, sequentially taking one 
Sample value at a time and carrying it through to increment 
the final integrator storage. Because of the limitation of 
only 4K of storage in the 160, the input data had to be buf- 
fered into the routine 120 values at a time and the output 
unloaded at 1500 values per run. This incremental method of 
printing output on the Cal-Comp plotter adds an extra line 
at each beginning point. 

Since the Hewlet-Packard sensor had a very smooth 
output compared to the Vibratron, it was used for all pre- 
liminary runs to test the tide cancellation of the filter 
section. Once smooth output was obtained from the filter by 
adjusting the precision and values of the weighting constants, 
the square and integrate routines were perfected on the 
Vibratron data. Whether the perturbations detected on the 
Vibratron tape were noise or actual events is of no concern 
to the detector's operation but is a problem of sensor 
design. 

Tracing the three interesting "events" detected on 
the Vibratron tape back to the original data, indicates 
pressure perturbations on the order of 4mm. Considering the 
tidal variations as the major noise source of the input data, 


and these 4mm. variations as signal, the overall improvement 
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in signal-to-noise ratio is on the order of 2.3-102 or 
67 db. A more detailed software analysis of system perfor- 


Mance, detection thresholds, power spectral densities, and 


evaluation 1s included in the Refs. 1 and 5. 
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VereeeNCLUSIONS 


The unusual event detector implemented with COS/MOS 
logic should provide a very useful and needed package to the 
growing number of tools used in oceanographic research. The 
detector can be used as an independent component in a larger 
system, or, due to its system structure design 1t can be 
used as the basis for a rather complex package. Several sub- 
systems for ascent and descent control, telemetry, data- 
storage, and sensor stabilization are actually needed to com- 
plete an underwater instrument capsule. The clock and timing 
circuits of the detector could be easily extended to cover 
other operations and simple gating of the memory package and 
addition of more memory cells, could make it a general buffer 


memory for the other functions attached to the system. 


ral 


APPENDIX A. COS/MOS 


LOGIC PACKAGES 


Lead numbers are not given since package styles vary. 


CD4001 QUAD 
TWO INPUT NOR GATES 





CD4011 QUAD 
TWO INPUT NAND GATES 


hZ 


CD4025 TRIPLE 
THREE INPUT NOR GATE 


CD4023 TRIPLE 
THREE INPUT NAND GATE 









SHIET REG@SIEk 





CD4015A DUAL 
FOUR STAGE REGISTER EIGHT BUFFERED INPUTS 


WITH SERIAL IN PARALLEL OUT 
CD4021A 8 STAGE SYNC 


PARALLEL INPUT/SERIAL OUT 













Ve Ve 
Vo V 
ENABLE G RESET 0 
VI VI 
Vpp ae Vo Vo 
4 STAGE COUNTER VT V7 
Vss a Vpp “SS 
a a 
CD4016 QUAD 
FIGHT DECODED OUTPUTS BILATERAL SWITCH 


Gp4022 piv iIpE sy EIGHT 
DIVIDER DECODER 


Us: 





eC 


7 STAGE 


REEPEE COunrtTeR 


So 





CD4009 HEX INVERTERS 
SEVEN BUFFERED OUTPUTS 


CD4009A 7 STAGE 
BUPPLE BINARY COUNTER 





CD4006 EIGHTEEN STAGE CD4013 DUAL D 
Soa Ce Shir Tl “RECS LER Peles! Lor 
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DATA INPUT . 


dresEr 


fai COMPLEMENTOR CIRCUIT USED SEVERAL TIMES IN DETECTOR 
DESIGN CAN BE FORMED FROM ONE CD4013 DUAL FLIP-FLOP AND 
ONE CD4011 QUAD TWO INPUT NAND GATE 


SWITCH 








ENS er 





LOS 2 OBES: 0 INPUT 


THE SERIAL SWITCH CAN BE FORMED FROM ONE CD4011 QUAD 
2 INPUT NAND GATE OR ONE-HALF OF A CD4016 QUAD BILATERAL SWITCH 
PLUS AN INVERTER 
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COMPUTER PROGRAM FOR SIMULATION 
OF THE UNUSUAL EVENT DETECTOR 


REWIND 1000 «+ MEMORY LOCATION AT BEGINNING 


OF FUNCTION BLOCK 
INITIAL 
OUT BUFFER / 1902 
ZERO 
IN BUFFER / 19°9 


LOAD 120 


PeleiGoa. iN 


1100 
emma COUNT CON- 
ie ie @ SECUTIVE ZEROS 
3 IF ZERO Ir 4 











FOUR#¥20 


D#0 
COMPLEMENT 
aS r 
TABLE 


=. a as 
ROUTINE . (ZERO ) 
ZERO 


700 . 740 


1550 : 7300 


IF PLOT ON INITIAL 
via eee 
MP T BUFFER 
QUT BUFFER es QU ) 
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7250 


1250 


7240 


0300 
0300 
0300 
0300 


nop 
nop 
nop 
nop 


idm 
jr 


ldi 
sbn 
Zjf 
jie 


1dn 
sti 


lad 
1s3 
jfi 


OPERATIONAL ROUTINES FOR VARIOUS ENTRY 
AND EXIT-OUTPUT OPTIQNS 


REMOVE ALL HALTS EXCEPT AT REWIND, OUTPUT, AND 


COMPUTATION COMPLETE 


TO READ OUT INPUT DATA 


1567 1600 
1ehe 0704 
12h5 1260 
1327 1330 
7250 2100 
o400 
POL. 
1560 


enter 
enter 
enter 
enter 


1dm 


a) Big E 


TO READ OUT FIRST DIFFERENCE 


1412 
1425 


0570 
T7200 
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enter 
enter 


1150 
0705 
1250 
7250 


truncated version 


0400 
Tako 





7200 
7300 
7400 


5 


to 


O1LOO 


GENERAL LISTING OF TSUNAMI DETECTOR PROGRAM 


TEMPORARY STORAGE LOCATIONS 
CONVERSION NUMBERS+ BINARY TO OCTAL 
SHIFT REGISTER STORAGE 

GENERAL INITIALIZE AREA 

ZERO INPUT ARRAY AND BEGIN TO READ 
READ 12ONUMBERS INTO 2000 to 3777 
BEGIN INPUT CONVERSION TO BINARY 
DIFFERENCE ROUTINE 

PREPARE FOR DOUBLE PRECISIGY 

CALL MULTIPLY SUBROUTINE 

SUM THE WEIGHTED DIFFERENCES 

SHIFT REGISTER IN 0400 to 0570 
SQUARE ROUTINE 

INTEGRATE ROUTINE 

INPUT BUFFER STORAGE 

OUTPUT BUFFER STORAGE 

MULTIPLY ROUTINE 

CLIPPER ROUTINE 

CAL-COMP PLOTTER ROUTINE 


RS-022 FLEX AND OCTAL READ ROUTINES 
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LISTING OF STORAGE LOCATIGNS SET ASSIDE FROM °°°° To 9°77 


into mit for 2 and 4 

temp storage of shifts 2 

ts of output of mt for 2 and 4 

ts value of p entering mult for 2 and 4 


ts of shifts 6 
ts 


ts for shift register routine 


O400 location of shift register words one 
O4OL 
0600 location of shift register words two 
0601 


aifference 1 
difference 31 
differrnce 61 
difference 91 
difference 121 


mult diff 1 
mult diff 31 
mult diff 61 
mmit diff 91 
mult diff 121 
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indigit + correction 

running sum first word 

running sum end word 

address of indigit being operated on 

kth digit of converted nos. and how many over 4 
7000 +~store 

count of nonzero and not 20 digits 

address of first of 4 indigits to be used 


reserve all remaining for rs-022 
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0000 


1000 


1010 


7101 
TOD 


jtL 


LOAD FLEX TAPE USING rso22 


REWUOYD TAPE AND BEGIN PROBLEM AT START 


initialize constants 


zero the flex read command for safety 
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1050 2200 
2000 


1060 


4100 
1056 
O400 
4100 


sbc 


njb 
ldc 


std 
hit 
gta 


ZERO INPUT ARRAY AND BEGIN TO READ 


halt when array is zero 
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1150 


1160 


1163 


1170 


ldc 
stm 
1dn 


std 
std 


BEGIN INPUT CONVERSIGN TO BINARY 


zero 61 and 62 


step 63 by one and read new digit 
if new digit is zero, step 5 and count 


zero the zero count 


test for 4 consecutive zeroes 


halt when out of input digits 


initialise constants 


go to zero erray and read more indigits 
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12350 


1240 


BEGIN TO CCNVERT TO BINARY 
if 20 step 4 


step 66 1 and go read another 


if it was 20 


are past 4 real nos 


zero 66 


read more 


zero 66 


store in 67 address of lst of 4 digits to be converted 
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1260 


1262 


13500 


1310 


1320 


BEGIN TABLE LOOK UP 


add new digit to running sum 


add 200 to indigit + corr. nos 


test to see if we have read 4 nos 


go back to 1260 and convert more 


separate nos into 2 words 


store 12 least simificant bits in shift 
register 0400 to 0570 


store most significant bits in 0600 to 0770 


halt input to binary is complete, now filter 
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1340 


1550 


1360 


1570 


1312 


sbi 


DIFFERENCE ROUTINE 


initialize constants in 20,22,24 


begin each difference here 


difference word 1, least significant digits 


for positive difference of word one 


for pos word 1 and pos word 2, zero the difference 


for pos word 1 and neg word 2, comlement the sign 


e 


for negative difference of word one 


for neg word 1 and neg word 2, zero the difference 
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1410 


1420 


for neg word 1 andpos word 2, complement the sign 


test for completion of all 5 differences 


step the difference ccunter by 36 


go back to 1340 and calculate another difference 


go to spike removal routine 
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1450 


1460 


1470 


1500 


1510 


std 
lde 


std 
lde 


std 
lde 


std 
ide 


std 
pel 


PREPARE FOR DOUBLE PRECISIQN 


for negative nos. store ones in 12 most significant 
bits 


for positive nos. zero most significant bits 


store miltiplication routine access location. 


load constant -3.919921875 


load constant 45.64179700 


load constant +1.000000000 


jump to call ARITH routine 


88 





1520 0101 
OO4 


1530 


1540 


pta 
Jpi 


CALL MULTIPLY SUBROUTINE 


miltiply difference 1 


miuitiply difference 31 


multiply difference 61 


mltiply difference 91 


multiply difference 121 


jump to swm the weighted differences 
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1550 


1560 


1570 


2047 
3051 
35055 
3055 
3057 
0300 
0300 
0300 
4100 
L000 
2101 
3600 
6770 
6203 
7101 
1600 
T7100 
7101 
73500 


SUM THE WEIGHTED DIFFERENCES 


store the result in 4000 


jump to the shift register 


halt when 4000 to 6770 is full of output 
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1600 


1610 


1620 


1630 


SHIFT REGISTER 0400 tO 0571 


shift data to right one sample in 0400 to 0570 


shift other half in 0600 to O770 


loop around again 


go back to convert another input nos to binary 


oul 





1555 hOSO std ENTER THE SQUARE ROUTINE 
7101 jri 


1700 2200 lde square routine 


2200 lde scale factor.....(O4#00 divides by 2) 


1710 0400 ldn zero most significant digits 


bo51 std 
4053 std 
4055 std 
O57 std 
2450 led rectify 
6302 njf 
4050 
al 
1720 0101 pta enter ARITH subroutine 
7OO# jp 
0031 multiply squares input 
0050 
0050 
0052 
0031 multiply scales output 
0052 
1730 0054 
0056 
2057 ldd load into 56 the squared and shifted results 
7101 By aut 
1560 
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1740 


1750 


1760 


1770 


7512 


7365 


LOOT 


ldc 


stm 


INTEGRATION ROUTINE 


input first value 


divide squared input by 8 


takes care of negative zero 


store the integrator output in the out. buffer memory 


increment memory one 


jump back into program 


corrections to CALCOMP PLOTTER ROUTINE 
for scaling 


corerections to initial section 
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7500 2200 lde CAL-COMP PLOTTER ROUTINE 


7510 2100 idm 


ooo 
koll std 
7101 spe 
7320 


7320 2011 ldd ccmpute difference from last point 


5410 sba 
4012 std 
6203 pjf 
7101 pig 
7340 

7101 glist 
75350 


7550 


7330 7402 write plus x one 


2012 ldd 
O701 sbn 
4o12 std 
6604 pjb 
7401 -X 
7101 Jf 
1350 


7340 7hOl write minus x one 


54.12 aod 
6702 njb 
7101 jPi 
7350 
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7350 0300 
7Tho4 


7360 


1370 


write minus y 


apes 


step input location one 


halt, output buffer is empty 


initialize storage number 4000. 


a5 





0101 


o141 


0161 


conversion for lst 3 octal digits of 4th indigit 


conversion for 3rd indigit 


conversion for 2nd indigit 


conversion for lst indigit 
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0501 


0321 


0001 
0003 


conversion forlast 2 octal digits of 4th indigit 


conversion for last 2 octal digits of 3rd indigit 
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ARITH 
A°closed subroutine for control of the CDC 168 Arithmetic Unit. 


PURPOSE ; 

This generator routine provides the capability to implement 22-bit 
fixed point arithmetic sequences using a three address code. Both 
fractional and integer arithmetic are provided. An indirect addressing 
option is included. 


USAGE: 

LOAD the ARITH subroutine anywhere in memory. (Space required 161¢ ) 
Store the load address in any convenient low-core directory cell (OOXX)> 
CALLING SEQUENCE IS3 


co ° 4 


PTA 
JPI XX OP Codes 

OP O1 = add 

A O02 - subtract 

B 30 = miltiply fractional 

C 531 - miltiply integer 

OP ho + divide fractional 

A hi - divide integer 

B 

C A and B: Operands 

° Cs: Result 
{NXT ) 

NOLES: 1. Any number of operations may be called for in sequence 


prior to return to normal 160 coding. 

2. Operands mist be in standard format for 22-bit fixed- 
point arithmetic (17s complement with allowance for 
Sign bit in each word). Least significant half-word 
in EVEN numbered cell, Most significant half-word in 
ODD numbrred mate. 


3. FOR normal direct addressing operands must be referenced 
only by even numbered address. 


4. Use of odd addresses will provide indirect addressing. 
5. Cells 70 - 77 are used as temporaries by the program. 


6. Timing - Approximately 900 microsecs. per operation. 
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ARITH 


SUBROUTINE FOR DOUBLE PRECISION 


OPERATICNS USING THE CDC 168 


init 


loop 


incr 


scode 
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sr-ARITH . 
control routine for 168 wit 
mic - 8/30/64 
02 
TT form address of lst argument 
70 set to Pill OP-A-B-C in 
76 10-71=-72=-T5 
oh 
75 
Th transfer arguments 


76 modify as required 
if odd(indirect) 


set to generate sel code 
33235 for plus 
3363 for minus 

70 3503 for mult 

O1 3502 for divide 





dvd 


mil) 


sub 


setout 


dvi 


hop 
dvf 


100 


set output addresses A-B 


these addresses OX for 
plus, minus, multiply 


process divide integer 


HOP links REPEAT to INIT 
process divide fractional 


common for divide 





t5 


repeat 
£45 


101 


sense sign of dividend 


fill neg zero for DD 


fill pos zero for DD 


activate 


skip unless divide 


set input address C 


skip unless MUF 


sense valid op code 
no valid - repeat 
valid - exit 
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